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Diabetes is one of the most common and critical diseases around the world, 
which need insulin injections to control the body’s glucose rate. A robust 
back stepping (BS) controller design based on fuzzy system is introduced in 
this paper to control the glucose level with the presence of meal disturbance. 
The controller’s design is based on Bergman’s mathematical model. 


Different fuzzy controller structures are implemented (fuzzy PI, fuzzy PD, 


and fuzzy PID) controllers along with BS controller named as (BS-fuzzy PI, 
BS-fuzzy PD, BS-fuzzy PID) controllers. Simulation results using 
MATLAB/Simulink show efficiency and robustness of the proposed design 
in terms of controlling the insulin concentration level in blood under meal 
disturbance and retaining the glucose level to its Basal value. 
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1. INTRODUCTION 

Diabetes is an illness caused by an insulin deficiency of blood. This causes an abnormal increase in 
blood sugar levels for a long time while the factors causing the disease are either hereditary or environmental 
[1]. Blood glucose is normally re gulated by the insulin hormone, which is produced by beta cells in the 
pancreas. In type 1 diabetes mellitus (T1DM), the beta cells are incapable of providing enough insulin due to 
the body’s attack (autoimmune destruction) and the glucose level becomes uncontrollable. This type needs to 
inject insulin into the human body to maintain a normal glucose concentration [2]. Hyperglycemia (Glucose 
concentration in bloodstream is higher than the normal level) leads to life-threatening of a patient. Unless the 
patient with diabetes mellitus is given sufficient insulin, his life will be at risk and makes the body vulnerable 
to many serious problems such as renal failure, nerve damage and limb amputation [3], [4]. A normal 
human's blood glucose concentration is in a narrow range (70-110) mg/dL. Diabetes is diagnosed when the 
human body is unable to control the normal interaction of glucose and insulin [5]. Therefore, insulin must be 
injected to regulate the level of glucose in the blood [6]. The closed-loop glucose regulation system is 
generally formed up of three main elements, which are: control method, insulin pump, and glucose sensor to 
estimate the insulin dose depending on glucose monitoring [7]. Figure 1 illustrates the block diagram of the 
insulin regulation system. Controllers are needed to set insulin treatment according to the measurements of 
glucose concentrations as shown in Figure 1. 

Various kinds of control units were introduced by many researchers in the literature. To name few, a 
second-order sliding mode control design is introduced to control the glucose concetration level [8]. Two 
robust control methods are presented in [9], to regulate the glucose-insulin system for Type I, LQ and Hoo 
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control are proved to be effective for providing acceptable control performance (By refusing to disturbances 
food meal). A new back-stepping based sliding mode control (B-SMC) approach is presented as a way to 
ensure practical glucose concentration tracking [10]. A back stepping controller design is proposed, based on 
Lyapunov stability theorm for the glucose-insulin model of type I diabetes in [11]. A comparison between 
back stepping sliding mode controller and the conventional state feedback controller is presented in [12]. A 
single model predictive controller to regulate insulin, and glucagon, is presented in [13]. An integrated of 
adaptive backstepping is proposed in [14], based non-linear controller with artificial pancreas, to control type 
1 diabetes. Three nonlinear controllers are proposed in [15]. Integral backstepping controller, backstepping 
controller and fuzzy logic controller, for the automatic regulation of the blood glucose level in diabetic 
patients. Under parametric uncertainties. A nonlinear explicit model predictive control, to regulate blood 
glucose in type-1 diabetic patient, is presented in [16]. A simulation model of glucose metabolism process as 
well as a backstepping linear quadratic Gaussian controller design, is proposed in [17], to regulate the blood 
glucose level in type-1 diabetic patients. 


Desired Blood Pump speed Insulin flow Actual Blood 


Glucose 5 Glucose 
—— 


Measured glucose 
concentration 


Sensor 


Figure 1. Closed loop insulin regulation system block diagram [8] 


In this paper, a fuzzy-PID, fuzzy-PI and fuzzy-PD controllers are combined with a back stepping 
(BS) controller and implemented to control the glucose level in the human body under meal disturbance. The 
paper is organized as follows into five sections. In section 2, the model is presented; BS design is explained 
and fuzzy logic system (FLS) is presented in section 3. Section 4 shows the proposed controller's efficiency 
with simulation results. Finally, a conclusion is presented in section 5. 


2. MATHEMATICAL MODEL OF BERGMAN'S SYSTEM 

Many different mathematical models of diabetic systems have been studied in previous works. The 
mathematical model presented in [18] is adopted in many researches considering the dynamics of glucose 
and insulin interaction in the blood system [18]—[20]. In this work, the mathematical model in [18] is used. 
The nonlinear differential equation is given by [20]. Physical meanings of the variables in (1) are shown in 
Table 1. 


G(t) = —piG(t) — (GCE) + Gg) X(t) + h(t) 


X(t) = -pX (Œ) + pY (t) 
Y(t) = —pa(V(t) + Ys) + i(t) /V, (1) 


Table 1. Physical meanings of the variables in Bergman's system equations 


Symbol Definition Unit 

G(t) Plasma glucose deviation mg/dL 
X(t) Remote compartment insulin utilization l/min 
i(t) Exogenous glucose infusion (control variable) mU/min 
Y(t) Plasma insulin deviation mU/dL 
h(t) The disturbance mg/dL.min 
Gp Basal glucose level 110 mg/dL 

Yg Basal insulin level 1.5 mg/dL 

Vi The insulin distribution volume 120 dL 
P1P2,P3,P4 The model parameters l/min, I/nin,dL/mUmin’?, l/min 


g(t) =x, > g(t) = X(t) 
x(t) = x2 > x(t) = X(t) 
y(t) = x3 > ý(t) = %3 (0) 
i(t) =u (2) 
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In (1), can be written as: 


x, (t) = —pıxı — (x1 + Gg)xX2 + A(t) 
X(t) = —P2X2 + p3X3 
X3(t) = —p4 (x3 + yg) + 1/V_u (3) 


the linearity form in the state space of the model was taken at the points of equilibrium with following 
specified values: 


Go = Xo = Yo =~Oh=0 (4) 
and for ig = p4YpV;, = 16.667. 


It is assumed that the state variable X(t) is also called slow variable, that is X(t): 


GCE) = —p,G(t) — (GCE) + Gp) X(t) + ACE) 
0 = -pX (t) + pY (t) 
Y(t) = —pa(V(t) + Ys) +i()/V, (5) 


xX) = EVO) (6) 
p2 
after applying in (6), the (1) becomes: 


G(t) = —p,G(t) — (G(t) + Gg) y(t) + h(t) 
Y(t) = =p4 (Y (t) + Ye) + i(t) /V, (7) 
This is due to the fact that the non-measurable X(t) variable (the use of isolated insulin) is a slow 


variable, which can be considered zero, thus eliminating the second equation of model (1) [21]. As a result, 
the linear low model is [22]. 


BSE ig p-p — 5, Go + | a + B i(t) + A h(t) (8) 
Since and, 

yi Gy =O (9) 

G(t) = x, Y(t) = x, (10) 


Le ee Tal de + lee 


1 0 0 0 
the system can be written as: 


iy = Pity — 7$ Gex +h (11) 


Xq = —p4X2 + 1/Viu (12) 


3. BACKSTEPPING CONTROLLER (BSC) DESIGN FOR BLOOD GLUCOSE 

The Backstepping approach is an effective tool for high-order system control problems [23], [24]. 
The goal of the proposed BS is (to stabilize the system when starting from any initial state and regulate the 
state to the origin (until reaching the origin). Let us consider the (11) as the first subsystem and the (12) as 
the second subsystem. Accordingly, consider in (11) with x1 as the state variable which is required to 
stabilize it via x2. x2 is considered here as a virtual controller, namely, we rewrite the (11) as follows: 


X, =V (13) 
ži = -p1 xı — = Gey +h (14) 
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in (14) is a first-order system with v as a control input. Let v be chosen as: 


v= Ben [-p1 xı + 4x] > v= Bes [—p, + A)lx, (15) 
where A > 0. Hence in (14) became: 

žı = —Ax, (16) 
the second state is: 

y =x Y (17) 

yH=x2-v 

ý = —paX2 + 1/Viu-—v (18) 


where: 


t = ay, i te 


P2 


according to the previous equations: 


A, = (a, ” + A) 


P2 


where the derivative v function is: 


Vahey 
v = Ay (Pid — CG axe +h 
Let u = V,[v — ay + p4x2] 


then in (18) becomes: 
y=-ay,a>0 


the output y goes exponentially asymptotically to zero as to. As y approach zero value x, > v. Finally, 
the control law is given by: 


u = ù — ay + p4x2 
u = V, |a (pix; — 2 Gaxa + h) — a(x, — v) + paxo] 


=V, [a Cp = rè GBX2 +h) — ax, — ad,x, + FA (19) 


where: 


ki =V, ri — ahı) 
k, =V, (AG) Ge — a +p,) 


by suitable selection for the parameters of the BS controller, these parameters are: 


à = 0.114 

A, = 0.15, æ = 0.02 
kı = —0.864 

ky = —1.536 


3.1. BS Controller based fuzzy system design 

BS controller based on fuzzy logic system (FLS) is widely used to solve the real-world problems 
due to its ability to compensate for structured and unstructured uncertainty. In this paper, a fuzzy PI, PD and 
PID are combined with BS controller to improve its performance. The proposed BS controller control law 
based on the fuzzy system is given by: 
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u(t) = —kıxı — k2x2 + Ufuzzy (20) 


where Uuzzy is fuzzy PI, PD, PID controller with: 


Up] = KpXy + ki f x2 (t) (21) 
d 
Upp = yxy + ka ZO (22) 
Uprp = KyXy(t) + ki f X2(t) + (kyz (©) + kg ZO) (23) 
PID ~ “p*1 i 2 prt d d(t) 


where k,,kaand k; are the proportional gain, the derivative gain, and the integral gain respectively. In this 
paper the standard procedure of fuzzy controller's design is followed. Generally, the structure of a fuzzy logic 
controller consists of three basic steps fuzzification, control rule base establishment, and defuzzification. 

In the fuzzification step, the physical values of the input and output are transformed into fuzzy 
linguistic sets. The fuzzy sets are normalized and arranged in subsets (intervals) according to the range of the 
input or output associated with membership function. These membership functions describe the degrees of 
the confidence of the input or output related to this range. This step is done to make the input physical values 
suitable with the fuzzy control rule base [25]. The membership functions considered in the design of the 
fuzzy PD, PI and PID controllers, are the same and are chosen as follows: i) fuzzy set 'plasma glucose’ has 
two memberships P (positive) and N (negative), ii) fuzzy set ‘plasma insulin' has two memberships P 
(positive and N (negative), and iii) the output controller is represented by three membership functions. And 
iv: Fuzzy set 'out' of PI, PD has three memberships P (out_positive), N (out_negative), and Z (out_zero) as 
shown in Figure 2(c). Fuzzy membership functions as shown in Figure 2. The Fuzzy memberships of Input 
for Plasma Glucose (g), Input for plasma insulin (y), and output (out) are illsetrated in Figure 2(a), Figure 
2(b), and Figure 2(c) respectively. 


-10 -8 -6 -4 -2 0 2 4 6 8 10 


Figure 2. Fuzzy membership functions: (a) fuzzy membership Input for plasma glucose (g), (b) fuzzy 
membership input for plasma insulin (y), and (c) fuzzy membership output (out) 


Fuzzy inference rules (FIR) are set of IF-THEN rules that control the relation between input and 
output of the fuzzy system and it is stated as follows: IF (set conditions satisfied) Then (deductive result set). 
These rules consist of a linguistic group associated with the input, output and control variables. It may also 
include fuzzy operators such as OR and AND [25]. In this work, the fuzzy rule systems for PI, PD and PID 
controllers are the same and it is represented by four rules as follows: 

- If g=N AND y=N THEN OUTEN. 
- If g=P AND y=P THEN OUT=P. 
- If g=P AND y=N THEN OUT=Z. 
- If g=N AND y=P THEN OUT=Z. 
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In the defuzzification step the fuzzy controller outputs are converted into crisp real signal (fou) and 
then send it as a control action to the physical plant [20]. There are several methods for defuzzification, the 
center of mass (COM) formula is used in this paper. 

_ LU(1)*output(1)+pu(2)*output(2)+......+uU(Nn)xoutput(m) 


pus pupu) + tuU) (24) 


Where: 
The control output as a crisp value is four. 
The number of rules is m = 4. 
The membership value of the output n is wu(m) 


4. SIMULATION RESULTS 

The proposed controller has been tested using MATLAB-Simulink (2016) software and its 
performance is assessed. The Simulink module is shown in Figure 3. The experimental results are carried out 
in a personal Computer where Windows 10 Pro is installed with 2.6 GHz processor and 8 GB of installed 
memory (RAM). 


Interpreted 

MATLAB Fen 

4 Interpreted | 
MATLAB Fen 


oH 


Table 


Figure 3. Simulink model of the proposed controller 


The controller has been applied to the diabetic model which is represented in (7) and the response of 
the system is tested with disturbance meal existence. Typical values considered for the parameters used in 
Bergman’s system equations are shown in Table 2. The parameters considered in the experimental results for 
the BS fuzzy controller are shown in Table 3. The behavior function for meal disturbance representing the 
exogenous glucose infusion is shown in Figure 4. 


Table 2. Parameters’ values considered in Bergman's system equations 


Symbol Definition Value 
Gp Basal glucose level 110 mg/dL 
Yp Basal insulin level 1.5 mg/dL 
Vi The insulin distribution volume 120 dL 
p1=0.028 I/min 
PhP p2=0.025 1/min 


wa The model parameters p3=0.00013 dL/m Umin? 


p4=5/54 l/min 
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Table 3. The parameters of proposed BS controller and PID fuzzy 
Parameter of BS controller 1 Parameter of BS fuzzy (PI, PD, and PID) 


A, = 1.5,a = 0.02 k,=1.4 
kı=-0.864 ka=0.5 
k=-1.536 ki=1.2 


be 


h(t) (mg/dl.min) 
ES 


0’ Å 
0 100 200 300 400 500 600 


Time(min) 


Figure 4. Disturbance meal function (exogenous glucose infusion) [9] 


Figure 5 shows two plasma glucose deviations. Figure 5(a) illustrates the simulation results for the 
plasma glucose deviation G(t) when BS-fuzzy controllers are used. The figures show the glucose deviation 
G(t) values does not exceed (120.45 mg/dl) when applying the BS-fuzzy PD and BS-fuzzy PID controllers 
which means a (9.5%) change over the basal value (110 mg/dl) while for a BS-fuzzy PI controller the 
maximum glucose level is around (123 mg/dl) i.e. (11.8 %) change over the basal value. 

Figure 5(b) shows a comparison between the simulation results for the plasma glucose deviation 
G(t) when BS controller is used and when BS-fuzzy PD is used. It’s obvious that when using the BS 
controller, the maximum glucose level reaches higher values of glucose deviation and hence permits higher 
percentage change over the basal value. Figure 6 shows the insulin deviation Y(t) when applying BS-fuzzy 
(PD, PI, PID) controllers. It is obvious from the figure that the insulin deviation does not rise over the Basal 
value for the insulin level (1.5 mU/dl). When applying the controllers this level is regulated to zero after less 
than 400 minutes. 

Figure 7 shows two insulin infusion controls. Figures 7(a) and (b) illustrates the insulin infusion 
control /(t) when applying BS, BS-fuzzy (PI, PD, PID) controllers. The insulin infusion is done to maintain 
the blood glucose level in the patient’s blood after having a meal. It is clear that the convergence time is less 
than 400 minutes. Higher infusion level can be seen for the BS-fuzzy PD and BS-fuzzy PID controller. 


14 120.4404 


PA ~ BS-fuzzy(pid) , ~~ BS-fuzzy(pid) 
122 -e—a —BS-fuzzy(pi) 120.4403 _ — —BS-fuzzy(pi) 
a >BS-fuzzy(pd) 
120 FN BSfuzzy(pd) sun ne 
3 f d 2 
BII8 \ E 120.4401 
ois = ma 
0 0 
ms ET] 
i 7 
" 1204398 
110 
0 100 200 300 400 500 600 131.998 132 132.002 132.004 132.006 = 132.008 132.01 
Time (min) Time (min) 
(a) 
124 
—Bs 
122 - - BS-Fuzzy (PD) 


0 100 200 300 400 500 600 
Time(min) 


(b) 


Figure 5. Plasma glucose deviations: (a) when applying BS-fuzzy (PI, PD, PID) and (b) when applying BS 
and BS-fuzzy (PD) controllers 
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Figure 6. Plasma insulin deviation when applying BS and BS-fuzzy (PI, PD, PID) controllers 
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Figure 7. Insulin infusion control: (a) when applying BS-fuzzy (PI, PD, PID) controller and (b) when 
applying BS only and BS-fuzzy PD controller 


CONCLUSION 
In this work, an improved BS controller based on different fuzzy controller structure was suggested 


to the nonlinear three-state minimal blood glucose model of Bergman, to control the blood glucose level. 
Simulation results are demonstrated through the Matlab tool, the efficiency of the proposed (BS and BS with 
fuzzy controller) under meal disturbance has obtained. BS and BS-fuzzy PI controllers show low 
performance as the glucose level exceeded the value of 123 mg/dL with lower insulin infusion as compared 
to the BS-fuzzy PD and BS-fuzzy PID. Both BS-fuzzy PD and BS-fuzzy PID show better performance as the 
maximum glucose level is less than 120. 5 mg/dL and higher insulin infusion. 
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